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ABSTRACT 

The velocity of the inner ejecta of stripped-envelope core-collapse supernovae (CC- 
SNe) is studied by means of an analysis of their nebular spectra. Stripped-envelope 
CC-SNe are the result of the explosion of bare cores of massive stars 8 M©), 
and their late-time spectra are typically dominated by a strong [O i] AA6300, 6363 
emission line produced by the innermost, slow-moving ejecta which are not visible at 
earlier times as they are located below the photosphere. A characteristic velocity of 
the inner ejecta is obtained for a sample of 56 stripped-envelope CC-SNe of different 
spectral types (lib, lb, Ic) using direct measurements of the line width as well as 
spectral fitting. For most SNe, this value shows a small scatter around 4500 km s _1 . 
Observations (< 100 days) of stripped-envelope CC-SNe have revealed a subclass 
of very energetic SNe, termed broad-lined SNe (BL-SNe) or hypernovae, which are 
characterised by broad absorption lines in the early-time spectra, indicative of outer 
ejecta moving at very high velocity (v ^ 0.1c). SNe identified as BL in the early 
phase show large variations of core velocities at late phases, with some having much 
higher and some having similar velocities with respect to regular CC-SNe. This might 
indicate asphericity of the inner ejecta of BL-SNe, a possibility we investigate using 
synthetic three-dimensional nebular spectra. 

Key words: 



1 INTRODUCTION 

Massive stars (> 8 Mg) collapse when the nuclear fuel in 
their central regions is consumed, producing a core-collapse 
supernova (CC-SN) and forming a black hole or a neutron 
star. CC-SNe with a H-rich spectrum are classified as Type 
II. If the envelope was stripped to some degree prior to the 



explosion, the SNe are classified as Type lib (strong He lines, 
and weak but clear H), Type lb (strong He lines but no H), 
and Type Ic (no He or H lines) ( |Filippenko|1997 |. 

Some CC-SNe, called broad-lined SNe (BL-SNe), ex- 
hibit very broad absorption lines in their early phase, result- 
ing from the presence of sufficiently massive ejecta expand- 
ing at high velocities. BL-SNe seem to be preferentially of 
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Type Ic [two exceptions are the Type lib SN 2003bg ( Hamuy 
et al.|2009|> and the Type lb SN 2008D ( |Mazzali et al.|2008| 



Modjaz et aL|2009| )]; se e also SN 1987K (jFilippenko||i988"] ), 
which is mentioned by ( Hamuy et al.|[2009 l. Some BL-SNe 
can reach kinetic energies of ^ 10 Oyl ergs. They are sometimes 
called hypernovae, and can be associated with long-duration 
gamma-ray bursts (GRBs) (see |Woosley fc Bloom|2006"l and 
references therein). However not all BL-SNe are associated 
with GRBs. 

An important question in the context of CC-SNe is how 
the gravitational energy is converted into outward motion 
of the ejecta during the collapse; see ( Janka et al.|2007 1 for 
a recent review. In GRB scenarios a relativistic outflow is 
launched by the central engine and deposits some fraction of 
its energy into the SN ejecta, probably preferentially along 
the polar axis, which might cause strong asymmetries (e.g., 
Maeda et al.|2002[ ). The ne arest, best-studied G RB-SNe are 
SN 1998bw / G RB 980425 flGalama et al.|1999[ ), SN 2003dh 
/ GRB 030329 |Matheson|2004| ), SN 20031w / GRB 031203 
( Malesani et al.|2004| , and SN 2006aj / GRB/XRF 060218 
( Pian et al.||2006 l, although it is not fully established that 
the GRBs (or X-ray flashes) accompanying nearby CC-SNe 
can be compared directly to high-redshift GRBs. CC-SNe 
may be characterised by asphericities although a jet does 
not necessarily for m jBlondin et al. 2003| Burrows et al 



20071 |Kotake et al.||2004| |Moiseenko et al.|[2006| |Takiwaki 
et al.|2009| >~ 

Extremely massive stars (> 100 Mq) are thought to 
end their lives as pair-instability SNe (PI SNe) . A star with 
enough mass to form a He core with more than 40 Mq will 
suffer electron-positron pair instability, leading to rapid col- 
lapse. This triggers explosive oxygen burning, leading to the 
complete disruption of the star ( Barkat et al.|196~ Heger fc 
Woosley||2005| >. This process can produce large amounts of 
st) Ni. The ejecta mass and the explosion kinetic energy are 
high, but the ejecta velocities are moderate (Scannapieco 
et al.|2005| >. 



Stripped CC-SNe, which lost part of their envelope be- 
fore collapse, offer a clearer view of their inner ejecta than 
SNe which have retained it. Thus, in this paper we exclu- 
sively address stripped CC-SNe (Types lib, lb, Ic); we do 
not include Type II SNe, despite using the term "CC-SN." 



Asphericities in their inner and outer ejecta (e.g., Maz- 
zali et al.|[2001 l are evident in at least some CC-SNc. Two 
main indicators arc velocity differences of Fc and lighter- 
element lines, and polarisation measurements (e.g 
19911. 



Hoflich 



Independent of their type, with time SNe become in- 
creasingly transparent to optical light, as the ejecta thin out. 
At late times (> 200 days after the explosion), the innermost 
layers of the SN can be observed. This epoch is called the 
nebular phase, because the spectrum turns from being ab- 
sorption dominated to emission, mostly in forbidden lines. 
In this phase the radiated energy of a SN is provided by the 
decay of radioactive 56 Co (which is produced by the earlier 
decay of 56 Ni). Decaying 56 Co emits 7-rays and positrons 
which are absorbed by the SN ejecta. As the deposition rate 
of 7-rays and positrons depends on the density and 56 Ni dis- 
tribution, the inner parts of the SN dominate the nebular 
spectra. 

Recently, several authors ( Maeda et al.|[2008| |Mo djaz 
et al.||2008 Taubenberger et al. 2009 1 have studied nebular 



spectra of CC-SNe. They concentrated on the shape of the 
[O 1] AA6300, 6364 doublet (which is produced by much of 
the mass) and concluded that torus-shaped oxygen distribu- 
tions might cause the double-peaked [O 1] profile observed 
in many CC-SNe nebular spectra. However, there is ongo- 
ing discussion regarding whether geometry is the dominant 



reason for this type of line profile (Mi lisavljevic et al.|2009 
also see Appendix A). 

Several authors have modelled nebular-phase spec tra of 
SNe to derive quantities such as the 56 Ni mass (e.g., [Maz- 
zali et al.||2004| |Stritzinger et al.|[2006| |Sauer et al.||2006 
Maeda et al. 2007a I, eje cta velocities (e.g., |Mazzali et al 
2007b|>, asp hericities (e.g., |Mazzali et al.|2005[|2007a||Maeda 



et al.J |2008[ ), and elemental abundances (e.g 



Maeda et al 

2007a|b| >. The nebular phase is especially suitable for study- 
ing the core of SNe. If different explosion scenarios are in- 
volved for different types of CC-SNe, one might expect the 
largest, most revealing differences to be in the central region 
of the explosion. Therefore, in contrast to the standard clas- 
sification of BL-SNe, which is based primarily on early-time 
spectroscopy and describes the velocity field of the outer SN 
layers, here we focus on the centre of the explosion. We have 
modelled the nebular spectra of over 50 SNe, the largest sam- 
ple of CC-SNe so far, and obtained a statistically significant 
representation of their core velocities. 

We describe our data set in Section 2 and the modelling 
procedure in Section 3. In Section 4 we test the reliability of 
the modelling approach. Results are discussed in Section 5. 



2 DATA SET 

We collected nebular spectra of 56 CC-SNe. This sample in- 



cludes all the spectra presented by Matheson et al. (2001), 



Modjaz et al. (20081, and Taubenberger et al. (2009) for 



which a spectral fit was possible. The most important crite- 
ria for selection were a reasonable signal-to-noise ratio and 
a spectral coverage of at least the region between 6000 and 
6500 A. Most spectra range from 4000 to 10000 A, allow- 
ing modelling of the Fe-group, oxygen, calcium, and carbon 
lines. If we found evidence for an underlying continuum, we 
tried to remove it using a linear fit. When several nebu- 
lar spectra were available for a given SN, we chose the one 
closest to 200 days, although the precise epoch has little 
influence as long as the spectrum is nebular, as shown in 
Section 4. 

Unfortunately, most of our spectra are not properly 
flux-calibrated. If an estimate of the 56 Ni mass was avail- 
able in the literature for SNe with uncalibrated spectra, we 
used these values (Table 1). For the remaining SNe with un- 
calibrated spectra we tried to estimate the 56 Ni mass from 
the light curve. We emphasise that the exact 56 Ni mass is 
not important for our study, as we show in Section 4. 

For several SNe only one light-curve point exists (gen- 
erally the detection magnitude in the V or B bands, or un- 
filtered). To obtain a (very crude) estimate of the 56 Ni mass 
from this single data point we employed the following pro- 
cedure. For the SNe with known 56 Ni mass and bolometric 
peak luminosity, the ratio of bolometric peak luminosity (in 
units of 10 42 ergs) to 56 Ni mass (in units of Mq) can be 
calculated. This ratio varies between 0.040 and 0.095, with 
a rather uniform distribution (see Table 1), which reflects 
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SN 


log(iBol) 


M m (M a ) 


P 


Mej.tot (M ) 


£ kin ,tot (W 51 ergs) 


References (-L pca k, Mm, -Mtot&Etot) 


19901 


42.4 


0.11 ± 0.02 


0.044 


3.7 ± 0.7 


1.1 ± 0.1 


1,1,1 


1993J 


42.2 


0.08 ± 0.02 


0.050 


~3 


~1 


2,2,13 


19941 


42.2 


0.065 ± 0.03 


0.041 


~0.9 


~ 1.0 


3,3,7 


199 /dq 


42.2 


U.15 ± U.Uo 


0.095 


~8— 10 


10—20 


5,5,5 


1997cf 


42.2 


n 1 q c; _l n not; 


0.085 


~8— 10 


10—20 


4,5,5 


199&DW 


42.8 


U.4y ± U.U4 


0.079 


~14 


~60 


3,3,15 


zUUzap 


42.1 


n (Tzq _L n no 
U.U/o ± U.Uz 


0.058 


2.5—5 


4—10 


3,3,14 


2003jd 


42.8 


0.3o ± 0.04 


0.057 


3 ± 0.5 


•7 + 3 


3,3,3 




42 4 


U.Z_L m U.UO 


n nsA 

U.UO 1 ! 


o.o o.u 


o.o y.u 


O, -LU, 1U 


zUUoaj 


42.7 


n on _L n n/t 
0.20 ± U.04 


0.040 




-^2 


3,12,12 


2007Y 


42.1 


0.06 ± 0.01 


0.048 


~0.5 


~0.1 


8,8,8 


2007gr 


42.2 


0.08 ± 0.02 


0.046 


1.5-3 


1.5-3 


10,11 


2007ru 


42.9 


0.4 


0.045 


1 3+ 1 ' 1 


c+4.7 
°-3.0 


9,9,9 


2008D 


42.2 


0.09 ± 0.02 


0.057 


~7 


~6 


4,4,4 


2008ax 


42.3 


0.08 ± 0.02 


0.040 


3-6 


~1 


6,6,6 



Table 1. Peak luminosities, 56 Ni masses, the ratio (3 = Mm (MqV-^BoI 42j the total ejecta mass, and the total kine tic energy for 15 
CC-SN e taken from the l iterat ur e (see referenc e s). Typ ic ally j3 w 0.06 with a di spersion of roughly 7 %. References: ^Elmhamdi et al. 
(|2004|) ^Bartunov et alTI (|1994| ) ^Valenti etaTl (|200Sa.|), ^Mazzali et al.| |2008) ^Mazzali et al.|l|2004| ), jPastorello et al.| (|2008[), jSauer 



et al. 12006) , IStritzinger et aL](|2009t , ^Sahu et al.|(|2009t, lu H unt er, in prep., 1] fValenti et al.|j2008b^ , ^Mazzali et al.| | |2006| ), ^Nomoto 
p993] . 'T 



et al 



^Mazzali et al.||2002[ l, ^Nakamura et al.| < p000^ , ^Taubenberger et al.| l |2006 



the variety of light-curve widths of different SNe (e.g., SNe 
with broader light curves have more 56 Ni for the same peak 
luminosity). From the estimates of the 56 Ni mass and the 
peak luminosity of these SNe, we can derive a relation 



M m (M Q ) 



-0.025 -^pcak,42- 



(1) 



This estimate should be compared with a similar one 
obtained for SNe la by Stritzinger et al.] ( |2006[ > , who found 
Mm (Mq) = 0.050 -Lpeak.42- The difference probably arises 
from the different densities and compositions of SNe la and 
CC-SNe. 

To estimate the peak luminosity from the measured 
magnitude, we first tried to determine the bolometric lumi- 
nosity at the time of detection. For some SNe an estimate for 
both the Milky Way and host-galaxy absorption is available 
in the literature. For most SNe, however, only the former is 
known ( [Schlegel et al.|1998[ ). In this case we assume a host- 
galaxy absorption between and 1.0 mag (unfiltered) and 
treat this range as an uncertainty affecting our estimate. 
If the Milky Way absorption is not known as well, we as- 
sume an uncertainty between and 1.5 mag (unfiltered). 
Unfiltered magnitudes are treated as bolometric, V-band 
magnitudes are converted to bolometric magnitudes using 
"iBoi = my — 0.3 ± 0.2 , and B-band magnitudes are con- 
verted using m,Boi = fna — 0.8 ± 0.4. For the distance mod- 
uli and the errors in the distance we took the values listed 
in NEE^for the SN host galaxies (Virgo+GA+Shapley). 
We then estimated the epoch of detection (which is close 
to maximum light for most of the SNe of our sample) and 
the uncertainty in this value from the references given in 
Table 2. Comparing to light curves of well-observed SNe we 
determined the value and uncertainty of the peak luminos- 
ity, including uncertainties related to absorption, conversion 



from filtered to bolometric luminosity, and the lack of a well- 
sampled light curve. Combining this estimate with Equation 
|TJ, we obtained 56 Ni masses for all 56 SNe of our sample. 
These are listed in Tables 1 and 2. The possible error of this 
method is very large, spanning roughly a factor of 20. This 
uncertainty estimate is very conservative; for most SNe the 
actual error should be much smaller. However, it is sufficient 
for our purposes (as we show in the Section 4). 



3 SPECTRAL MODELLING 

The spectra were modelled using the nebular code of |Maz-| 



zali et al. (2001 2007a I in the stratified version. A Monte 



http : I /nedwww.ipac. caltech.edu/ forms/byname. html 



Carlo routine is used to calculate the deposition of energy 
(which is carried by the 7-rays and positrons produced by 
the decay of 56 Ni and 56 Co) in each shell. The gas heat- 
ing caused by this process is balanced by cooling via line 
emission. The excitation and ionisation state of the gas, to- 
gether with the electron density and temperatures, are then 
iterated until the line emissivity in each shell balances the 
deposited energy. The emission rate in each line is calcu- 
lated solving a non-LTE matrix of rates to obtain the level 
populations | Axelrod|1980 1. 

The general principle of such modelling is that different 
velocity shells are characterised by different element abun- 
dances and densities, which lead to different ratios of line 
fluxes in each shell. The total emissivity of a shell is con- 
trolled by the energy deposition, which depends on the den- 
sity and the 56 Ni mass of the shell and the neighbouring 
ones. Therefore, each shell has a certain emissivity inte- 
grated over all wavelengths and a characteristic line pro- 
file caused by the Doppler shift. Each line in the emerging 
spectrum is the result of the superposition of single compo- 
nents from each shell, with different characteristic widths. 
A change of the emissivity of one line results in a variety of 
profile variations of other lines which have to be modelled 
iteratively until the complex shape of the full spectrum is 
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SN 
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Ay (mag) 


m — M (mag) 
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log(ipoak) 


M Ni 
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1983N 


11. 3y 
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1991L 
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13.9jv 
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5 ± 5 


42 53+ - 56 


20+ 1 06 


36 * 


1995bb 












0.2 


37 


1996D 


18. 2y 


0.509 
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42.10+?, 68 

u 0.14 


0.07 + n'n? 

w w 1 —0.04 


38,* 


1996N 
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29 


1996aq 


14.7y 
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32.20 


± 4 




0.24+i' 48 


7 * 


1997B 


16.5jv 


0.243 


33.17 


10 ± 3 


42.40+n « 

— 0.06 


0.15+1 
—0.07 


39,* 


1997X 


13.5jv 


0.091 


31.15 


4 ± 4 


42.61+S-S 

— 0.00 


0-25_q ^2 


40,* 


2000ew 


14.9at 
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30.16 


14 ± 7 


41.88+",!!" 

— U.Uu 


0.05+0'^ 

— 0.0.2 


41,* 


2001ig 












0.13+S-jg 


8 


2003bg 


15.0jv 
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31.24 


14 ± 7 


42.25+°'«° 

— U.Uu 


0.11 " 

— U.U7 


10,* 


2003dh 












0.1 


9 


2004ao 


14.9at 


0.348 


32.35 


± 4 


42.56 '5 4 nR 

— U.Uu 


0.22+ n ?? 

— U.ll 


13,* 


2004dk 


17.6jv 
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32.15 


-10 ± 10 


41.87+"'^ 

— 0.06 


0.044+0^ 8 

— O.Ozl 


15,* 




lo.ojv 


n in 
U. 1U 


oi.Uz ^virgoj 


U ± 




U.22_ n 




2004gq 
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-4 ± 3 
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14.9jv 
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—3 ± 5 
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2004gv 
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— O.Oo 


26 


2005kl 
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>> 1 
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2006F 


16.7jv 


0.629 


33.70 


7 ± 7 


42.63±g;~ 




19,* 


2006T 


17.4jv 


0.246 


32.68 


-11 ± 2 
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— U.14 


0.07+ ,' 383 

— U.U4 
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2006gi 
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3 ± 5 
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0.06 


0.11+0'^ 

—0.06 
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2006ld 
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0.057 


33.74 


9 ± 4 


42.73+°-" 


0-33+|«° 


18,* 


2007C 


15.9iv 


0.140 


32.15 


2 ± 4 


42.03±°$ 


0.065iU; 8 ^ 


20,* 
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18.0at 
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34.75 ± 0.25 


10 ± 6 
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— U. 1U 
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— 0.0 7 


21,* 


2007bi 


18.3iv 
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43 64+0-86 

1O.D*_ og 
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z -°-1.3 
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2007ce 


17A N 


0/1.5 
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43 12+ ' 72 




43,* 


2007rz 


16.9at 
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42 52 +(K(i() 
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41 84+0-60 
tJ-O*_ .20 
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44 
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Table 2. Observed magnitude (V-band, B-band, N = unfiltered CCD), extinction (in the V band), distance modulus, epoch, estimated 
peak luminosity, and 56 Ni mass for 41 CC-SNe, together with references for the magnitude and absorption. For SNe 1987M, 2001ig, 
2003dh, and 2005bf, estimates for the 56 Ni mass or luminosity are available and referenced. The uncertainty on the epoch was taken 
from |Modjaz et al,| l |2008| ) and [Taubenbcrgcr ct al. (20091, if available, or estimated from the detection report. SNe without estimate for 
absorption are referenced with "*." For SNe 2007ce, 20071, and 2007bi, we have no information about the host galaxy, so we assumed 
an absorption between and 1.5 mag (unfiltered). SN 2004gk is a member of the Virgo cluster and we ther efore took the V ir go distance 
modulus for this SN, as no reliable estimate for the distance of the host galaxy is available. Refere nces: ^Tsvetkovj l|l986|), ^Begelman 
Sarazin|l|l986|), ^Nomoto et al ||l990|l, jClocchiatti fc Wheeler| (l997|l, ^fComez fc Lopez] ( |l994t, ^Elias et a" jl990\ , iNakano et ah 
' ^Silverman et al.| (|2009aJ, ^Deng et al | (|2005[ ) , 1 HWood-Vasey fc Chassagne| ||2003| l 1J |Monajd et al.| 2004|, 1 1Maund et al. 
HSinger fc Li| |2004|l HQuimby et al.| j2004| , 1 E |Graham fc Li| ||2004|l, i lMonardTj2006L lv |Itagaki et al.j \2006\ , ^Frieman 
.20061, 19 pimai fc Villi] (|2006|, ^Puckett et "at] pOpTt , 2 jjm et al.| (|2007t, ^jDimai fc Migliardi|(2005t, ^^Nugentl <2007> , 2 ^Parisky fc, 
|Li| ( |2007| l, ^ ^Blondin fc Calkins| ( |2008[ l, ^Maeda et al~] ||2007b|), ^ jMonard et al.| | |2004| , 2t |Taubenberger et al.| j2005| ~but absorption 
h ighly uncertain, as sume a aB Ni mass of 0.2 M Q , at jWilUams et aL] l |1996| ) but no light-curve information, assume a ati N i mass of 0.2 Mq, 
3 lPugh et al-M2004t bu no distance information available, assume a a0 Ni mass of 0.2 M , 31 |Clocchiatti et al.|(|l996|l, 3 ^Perlmutter et al. 
(1990) , 3 ^McNaught et al.| ( |l99l[ l but no light-curve information available, assume a 56 Ni mass of 0.2 M , ^Clocchiatti et al.| ( |2001[ | 
but distance and reddening highl y uncertain, assume a 56 Ni mass of 0.2 M 3; ]PoIIas fc Maury| jl99l] but no light-curve information, 
assume a 56 Ni mass of 0.2 M , 3i ]Filippenko fc Korth| (|l99l|l, 37 |Tokarz et al.] jl995> but no light-curve information available, assume a 
56 Ni m ass of 0.2 M , 3 8 jDrissen et al.|(1996) , 3t jGabrijelcic et al.|^1997) , ^Nakano et al.|fl997] , 4 1 |Puckett et al.|(|2000t , 4 ^Schmidt et al.| 
|2005| but no light-curve information, assume a 5fi Ni mass of 0.2 M , 4 ^Quimby et al.| (2007) , 4 jBrown et al.|(2008) but no light-curve 
information and no detailed epoch information, assume a 56 Ni mass of 0.2 M . 
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SN 


Type 


Epoch (days) 


ID 


v a (km s 1 ) 


i>50 (km s 


1983N 


lb 


226 


Y 


3797 


2630 


1985F 


Ib/c 


280 


Y 


4920 


2456 


1987M 


Ic 


111 


Y 


5486 


3701 


1990B 


1c 


110 


N 


5405 


5091 


19901 


lb 


237 


Y 


4828 


2899 


1990U 


Ic 


181 


Y 


3488 


3021 


1990W 


Ib/c 


183 


Y 


4803 


3425 


1990aa 


Ic 


141 


? 


4216 


4368 


1990aj 


Ib/c 


150-250 


? 


5122 


3034 


1991A 


Ic 


177 


Y 


5262 


3636 


1991L 


Ib/c 


100-150 


N 


4140 


3316 


1991N 


Ic 


274 


Y 


4278 


3239 


1993J 


lib 


205 


Y 


4029 


3070 


19941 


Ic 


147 


Y 


5057 


3967 


1995bb 


Ib/c 


150-400 


Y 


5154 


4410 


1996D 


Ic 


214 


Y 


5228 


3624 


1996N 


lb 


224 


N 


3736 


3047 


1996aq 


lb 


226 


N 


5846 


3451 


1997B 


Ic 


262 


Y 


4801 


3317 


1997X 


Ic 


103 


Y 


4680 


3420 


1997dq 


BL-Ic 


217 


Y 


4594 


3361 


1997ef 


BL-Ic 


287 


Y 


4681 


2733 


1998bw 


BL-Ic 


201 


Y 


6340 


3602 


2000ew 


Ic 


122 


N 


4184 


3001 


200 lig 


lib 


256 


N 


5027 


3241 


2002ap 


BL-Ic 


185 


Y 


6219 


3729 


2003bg 


BL-IIb 


279 


N 


4736 


3205 


2003dh 


BL-Ic 


229 


? 


4342 


3085 


2003jd 


BL-Ic 


317 


N 


6850 


5593 


2004ao 


lb 


191 


N 


4555 


3158 


2004aw 


Ic 


236 


Y 


5007 


3234 


2004dk 


? 


333 


Y 


5465 


4338 


2004gk 


Ic 


225 


Y 


4623 


3171 


2004gq 


lb 


297 


Y 


6697 


3039 


2004gt 


Ic 


160 


N 


4513 


3371 


2004gv 


Ib/c 


299 


Y 


4792 


3158 


2005N 


Ib/c 


70-120 


N 


4001 


3591 


2005bf 


lb 


209 


N 


3864 


3628 


2005kl 


Ic 


160 


Y 


5074 


3365 


2006F 


lb 


314 


N 


4491 


2966 


2006T 


? 


371 


N 


4202 


4436 


2006aj 


BL-Ic 


204 


Y 


6540 


5100 


2006gi 


? 


148 


Y 


4589 


3180 


20061d 


lb 


280 


N 


4086 


3182 


2007C 


lb 


165 


N 


4787 


3520 


20071 


BL-Ic 


165 


N 


6085 


3277 


2007Y 


lb 


270 


Y 


4331 


3025 


2007bi 


PI? 


360 


Y 


5487 


3756 


2007ce 


BL-Ic 


310 


Y 


6172 


4461 


2007gr 


Ic 


158 


Y 


4480 


3228 


2007ru 


BL-Ic 


200 


Y 


5811 


3981 


2007rz 


Ic 


292 


? 


4998 


3785 


2007uy 


lb 


111 


Y 


6103 


4563 


2008D 


BL-Ib 


86 


N 


5847 


4040 


2008aq 


lib 


>130 


Y 


4119 


2885 


2008ax 


lib 


246 


N 


4100 


2821 



Ref. 



Silverman et al. (2009b) 



Hamuy et. al. (in prep.) 
Bersier et. al. (in prep.) 



Modjaz et al. (20081 



Modjaz et"al~| [2008 ) 
Modjaz et al. (2008) 
Modjaz etlSl pOOSl 



Modjaz et al. (20081 



Modjaz et al. 



Modjaz et al 
Mazzali 



(2008) 
I2U0HI 



Stritzinger et al. I 2009 1 
Mazzali 
Matheson 



Valenti et al. 


(2009 


Sahu et al. 







Mazzali 
Mazzali 
Matheson 
Navasardyan et. al. (in prep.) 



Table 3. SN type (classified by early-phase spectra), the epoch of our spectra (relative to maximum light), and the characteristic 
velocities v a and 1150 both corrected by 5% per 100 days. The maximum errors of v a and D50 are estimated to be ± 14% and ±10%, 
respectively. SNe marked with "Y" are fit by our ID shell modelling quite well. For SNe marked with "N" the central parts of the line 
profiles are not reproduced well and therefore probably could be improved with multi-dimensional modelling (v a would not change much, 
as v a is dominated by the outer parts of the line profiles, where most of the kinetic energy is located). For SNe marked with "?" the 
spectra arc too noisy to categorise as "Y" or "N." In total we have 12 BL-SNe (2 of them are GRB hypernovae), 17 regular SNe Ic 
(27 with BL-Ic), 12 regular SNe lb (13 with BL-Ib), 4 regular SN e lib (5 with BL-IIb), 1 po ssible PI SN, and 10 SNe for which the 
classification is unclear (they are certainly of the CC-SN type). See |Taubenberger et al.]||2009[) for spectra referenced by "*." 
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reproduced. As a complicated structure must be found in 
order to produce a certain spectrum, the method is quite 
reliable in determining the velocity field of the ejecta. 

We start modelling each spectrum using a CC-SN model 
used by |Mazzali et al.| ( |2002) for SN 2002ap. Such a model 
contains all of the information about density, mass, velocity, 
and element abundances. We correct the 56 Ni mass to the 
values listed in Tables 1 and 2, and we scale the synthetic 
spectrum to match the observed one. 

A clear modelling of the Fe-group lines is impossible 
for most nebular spectra because these lines are generally 
weak and are therefore affected by noise and background. 
The oxygen line is therefore taken as a tracer for the 56 Ni 
distribution. This means that our 56 Ni zone extends out to 
the point where the oxygen line can no longer be separated 
clearly from the background. Hence, a narrow oxygen line 
results in a more central 56 Ni distribution in our models, but 
it is important to note that this assumption may not exactly 
reflect the situation in real SNe. This may cause some epoch- 
dependent error which we try to quantify by comparing the 
time dependence of line widths in our models and the ob- 
served line widths (see Section 4). Above the 56 Ni/0 zone 
we set the density to zero, as this region cannot be probed 
with the nebular approach. We discuss this in more detail 
in Section 4. 

Other elements, such as O, Ca, and C, are distributed 
over the velocity shells until the line profiles are matched. 
We pay special attention to fitting the O line, since oxygen 
is typically the most abundant element in stripped-envelope 
CC-SNe. 

The modelling process yields the abundances, masses, 
and velocities which best reproduce the observed spectra 
(see Figures [l] and [J for two examples). The mass and veloc- 
ity distribution can then be integrated to obtain the total 
mass and kinetic energy of the model. The ratio of core ki- 
netic energy (in units of 10 50 ergs) to core mass (in units of 
Mq) is termed a in this paper: 



E, 



M cj (M e 



(2) 



This parameter is measured for all 56 CC-SNe of our sample 
and is converted to a characteristic velocity (see Table 3), 



2af0 4 



i/2 r 



[km s 



(3) 



The largest uncertainty in the estimate of a is caused 
by the background. The kinetic energy is dominated by the 
outer parts of the oxygen line and these are superposed on 
background lines and noise, which are difficult to distinguish. 
To quantify the uncertainty that this causes, we tried dif- 
ferent fits for several spectra (see Figure [l] and [5] for two 
examples). We found that, depending on the background 
assumed, variations of up to ±15% in a can occur, which 
translates into an uncertainty of ~ 7% in v a . An exact treat- 
ment of this problem is difficult, as background subtraction 
is arbitrary to some degree. Therefore, based on this uncer- 
tainty, we estimate an error of 7% in v a . 

As we show in Section 4, the line width slowly evolves 
with time, causing v a to decrease by roughly 5% every 
100 days. To handle this problem, we corrected v a by 



(1 + 0.05(i(days) — 200)/100) and assume that this causes 
an additional error of ±5% in v a . 

In addition to v a , we also measure the half width of 
the oxygen doublet [O i] AA6300, 6364 at half-maximum 
intensity (HWHM), A 50 , for all SNe of our sample. This 
can be directly converted to a velocity via the redshift 



«50 = 



A 50 
6300 ' 



(4) 



where c is the speed of light and A50 is given in A units. 
This is another characteristic velocity of the SN core ejecta, 
which can be compared to v a (see Figure |3j. There seems to 
be a linear trend between v$o and v a , as expected; exceptions 
are mostly caused by specific features of some line profiles, 
as discussed in detail in Section 4. 

Estimating v$o has the advantage that it can be done 
easily and it does not require modelling. The disadvantage 
is that u 5 o contains some small contribution of the [O 1] 
A6364 line which causes some error. Since the ratio of [O 1] 
A6300 to [O 1] A6364 is ~3 in the nebular phase, this error is 
small; it depends on the exact shape of the individual lines, 
but it should be <20% in the worst case, as we found by 
superposing Gaussians. More importantly, it is influenced 
by the shape of the inner line profile (which in some SNe 
has a double- or even triple-peaked shape) much more than 
by v a . In addition, V50 measures the velocity at one point 
(half height), while v a is the result of an integration over 
the entire line profile and includes the nonlinear weighting 
of different emission regions. The background causes some 
uncertainty in the estimate of V50. We measure V50 for several 
SNe while varying our assumptions about the background, 
and estimate an uncertainty of ~5%. Moreover, we apply 
the same time-dependent correction as for v a and assume an 
additional error as for v a (±5%). Adding the possible errors 
due to background, 56 Ni mass estimate, and time evolution 
of the line width, we estimate a maximum error of ±14% for 
v a and ±10% for D50, which does not necessarily mean that 
i>50 gives a more correct estimate of the core velocity. 



4 DISCUSSION 
4.1 Tests 

In order to ensure that the estimates of v a are reliable, it 
is useful to verify that the data do not show behaviours 
that are not included in our modelling. Since the spectra 
that we used were obtained at very different SN epochs, 
it is vital that v a does not depend strongly on SN epoch. 
Figure [4] shows the time evolution of v a for three SNe with 
high-quality nebular spectra at several different epochs. The 
temporal evolution of v a is in fact weak, at most ~10% 
over 200 days, which is comparable to the general modelling 
uncertainty. Figure [5] shows the corresponding evolution of 
the line profiles. In Figure [6j v a is plotted against SN epoch 
for the entire sample; no strong time dependence is seen. 
For V50 the situation is similar. However, as can be seen for 
example in SN 1998bw (see Figures [4] and [5| , small features 
in the line profile can have disturbing effects on vso- The 
influence of epoch on the characteristic velocities will be 
more fully discussed in Section 4.2. 

To test whether the often highly uncertain estimate of 



© 0000 RAS, MNRAS 000, 000-000 



Characteristic Velocities of Stripped- Envelope Core- Collapse Supernova Cores 7 




6000 



Figure 1. Nebular spectrum of SN 1990U (v a = 
3450(blue)/3633(red) km s" 1 ) at t = 184 days (black) (up- 
per panel) and the oxygen doublet [O i] AA6300, 6364 isolated 
(lower panel). Two models are shown in blue and red, illustrating 
the uncertainty in the background. This spectrum is the one 
with the lowest v a in our sample (narrow oxygen line). The 
characteristic velocity is v$o = 2996(blue)/3124(red) km s -1 . 
The horizontal lines show the full width at half-maximum 
intensity (FWHM). 




5000 6000 7000 

Wavelength [A] 



1 




6100 6200 6300 6400 6500 

Wavelength [A] 



Figure 2. Nebular spectrum of the BL-SN 2006aj (v a = 
6542(blue)/6870(red) km s" 1 ) at t = 204 days (black) (upper 
panel) and the oxygen doublet [O i] AA6300, 6364 isolated (lower 
panel). Two models are shown in blue and red, illustrating the un- 
certainty in the background. This spectrum has one of the highest 
v a in our sample (broad oxygen line). The characteristic velocity 
«50 = 4992(blue)/5190(red) km s . The horizontal lines show 
the FWHM. 

the 56 Ni mass influences v a , we alternatively increased and 
reduced the 56 Ni mass by a factor of five, therefore spanning 
a factor of 25 in 56 Ni mass for several randomly selected SNe 
of our sample. This corresponds roughly to the maximum 
uncertainty in the 56 Ni mass estimates. We found that v a 
depends only weakly on 56 Ni mass: the difference is always 
less than 2%. For example, when changing the 56 Ni mass of 
SN 1987M by a factor of 25, from 0.6 to 0.024 M©, v a in- 
creased by only 1.3%, which is considerably smaller than 



Figure 3. The characteristic velocities V50 [km s — 1 ] vs. (v a — 
1500 km s _1 ). Both characteristic velocities seem to correlate (as 
expected). The offset by ~1500 km s" 1 is caused by the different 
velocities cutoffs used to calculate U50 (half height) and v a (full 
line profile). BL-SNe are shown in red. SNe Ic are marked with 
blue, SNe lb with green, and SNe lib with red dots. Even taking 
into account the substantial errors, it is obvious that BL-SNe have 
a large scatter of characteristic core velocities. The data point at 
the upper right is SN 2003jd (see Table 3). The four BL-SNe 
at low characteristic velocities are SNe 1997ef, 1997dq, 2003dh, 
and 2003bg. The intermediate characteristic velocity BL-SNe are 
(from the left to the right) SNe 20071, 1998bw, 2002ap, 2008D, 
2007ce, and 2006aj. The three BL-SNe 20071, 1998bw, and 2002ap 
fall off the v a /v$Q relation owing to the strong broadening toward 
the base of their [O 1] AA6300, 6364 line. 

other modelling uncertainties. Figure [7] shows v a against 
56 Ni mass; there is clearly no correlation. We conclude that 
the maximum error introduced by the uncertainty in the 
56 Ni mass is ±2%. The independence of v a on 56 Ni mass is 
discussed in more detail in Section 4.2. 



4.2 Discussion of the Method 

We obtained characteristic velocities for 56 CC-SNe. Before 
we turn to the results we discuss our methodology. 

In Section 2 we tried to determine the 56 Ni mass for 
29 SNe in our sample without calibrated spectra. A typical 
uncertainty in 56 Ni mass is a factor of 5-10 up and a fac- 
tor of 2-3 down. The uncertainty upward is mainly caused 
by the uncertainty in absorption and epoch. There are also 
uncertainties in distance, magnitude, converting magnitude 
to bolometric luminosity and finally to 56 Ni mass, which 
are smaller. Our 56 Ni masses might be systematically over- 
or underestimated and one should be cautious when using 
these values. The 56 Ni masses listed in the literature (18 
SNe) have an average of 0.176 M , while the SNe with 56 Ni 
masses estimated in this paper have an average 56 Ni mass of 
0.225 (28 SNe, excluding 9 SNe where we could not estimate 
a 56 Ni mass and excluding a possible PI-SN). Considering 
the relatively small numbers, the agreement is reasonable 
but might indicate a small systematic over-estimate of 56 Ni 
mass. In any case, as shown in Section 4.1, an uncertainty 
of a factor of 25 in 56 Ni mass will not influence our results 
significantly. 
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Figure 4. The [O i] AA6300, 6364 line profiles (left to right) of SN 
2002ap (day 156, black; day 185, blue; day 237, violet; day 274, 
green; day 336, yellow; day 386, red), SN 1993J (day 205, black; 
day 236, blue; day 255, green; day 295, yellow; day 367, red), and 
SN 1998bw (day 201, black; day 337, red; day 376, blue). The 
temporal evolution is quantified in Figure [5] 
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Figure 6. The epoch (days after maximum light) vs. character- 
istic velocity v a (km s _1 ) for all SNe at epochs listed in Table 3. 
SNe Ic are shown in red, SNe lb in blue, and SNe lib in green. 
SNe without type classification are shown in black. The effect of 
temporal line-width evolution is weak. 
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Figure 5. Temporal evolution of v a (full line) and nso (dashed 
line) for SNe 2002ap (red), 1993J(green), 1998bw (blue), 2007uy 
(light blue), and 2007gr (black). One can see that variations are 
rather small, about 5% per 100 days. The rapid drop of D50 of 
SN 1998bw is caused by the "knot" at half height of the oxygen- 
line profile, which "drops below" half height in the late phase. 
The increase of line width with time in SN 1993J must be some 
background effect (note the rise of the red wing at late times, due 
to Ha emission) which could not be subtracted. For SNe 1993J, 
1998bw, and 2002ap, similar temporal evolution was found by 
|Taubenberger et al.| | |2009| . 



We checked for possible correlations between 56 Ni mass 
and Va - This could be important in two different ways. First, 
the uncertain estimate of the 56 Ni mass might introduce 
some error in the determination of v a . This is not the case, 
however, as shown in Section 4.1. Increasing the 56 Ni mass 
by a constant factor throughout the ejecta will increase the 
energy deposited in any shell by the same factor; hence, 
the relative contribution of each shell to the emission will 
remain constant. As the emitted energy per particle also 
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Figure 7. The logarithm of 56 Ni mass (in solar mass units) 
against the characteristic velocity v a . No relation is seen for our 
full sample. There might be a weak correlation for the SNe with 
56 Ni mass estimates taken from the literature (shown in red), 
where the errors in the 56 Ni mass are much smaller. The mean 
56 Ni mass of observed CC-SNe seems to lie around 0.20 M (0.176 
Mq for 56 Ni mass estimates taken from the literature, 0.225 Mg 
for the 56 Ni mass estimates in this paper), with a broad diversity. 
Neglecting the possible PI-SN 2007bi, 56 Ni masses range from 
0.05 Mq to 0.80 Mq . Because of the large uncertainties of our 
estimates, these upper and lower limits are inaccurate. 



remains rather constant (the mass of other elements must 
be increased accordingly to the 56 Ni mass), the spectrum 
of each shell does not change much. Therefore, the emerg- 
ing spectrum is nearly constant apart from small differences 
arising from the small shift of line ratios in each shell (result- 
ing, for example, from changes of the ionisation balance). 

Second, there could be a physical correlation between 
56 Ni mass and core ejecta velocities. For our full sample, this 
does not seem to be the shown in Figure[7] However, 

our large uncertainties in the 56 Ni mass make a stringent 
conclusion impossible. There might be a weak correlation 
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for the subgroup of 18 CC-SNe with Ni masses taken from 
the literature (the uncertainties in these estimates should be 
much smaller), but the situation is not entirely clear owing 
to the substantial scatter. 

Given that in most spectra we cannot determine the 
56 Ni velocity very accurately (in CC-SN nebular spectra Fe- 
group lines are usually weak, strongly overlapping, and can- 
not be easily separated from the background), we can draw 
no conclusion about possible relations between 56 Ni mass 
and Fe-group element velocities, which might differ from 
light-element velocities in aspherical SN models. However, 
there is no correlation between the total kinetic energy or 
the total ejecta mass and the characteristic core velocity (for 
the 15 SNe listed in Table 1). A weak correlation of the ra- 
tio (i? e j,tot/A/ e j,tot) 1//2 with the characteristic core velocity 
is found, albeit with large scatter (see Figure [8]). 

The epoch at which a SN reaches its nebular phase gen- 
erally depends on the mass and the ejecta velocity. Some of 
our spectra are quite early (~90 days after maximum light), 
though they all seem to be sufficiently nebular to be treated 
with our modelling. The epochs of our spectra vary between 
100 and 400 days (after maximum light), and we have shown 
in Section 4 that this large span of epochs will not influence 
our results much. From detailed modelling of five SNe with 
good spectral coverage in the nebular phase (extending over 
200 days), we can estimate that the characteristic velocity 
may decrease by ~5% every 100 days. This is caused by the 
decreasing importance of the outer layers, which become less 
luminous. 

In principle, if the behaviour of the positrons and the 
detailed distribution of 56 Ni were known, we could repro- 
duce this line-width evolution. However, since these param- 
eters are unknown, it is necessary to make some assump- 
tions about both properties. As described in Section 3, we 
assumed that both the 56 Ni distribution and the positron 
deposition trace the oxygen line profile, leading to an oxy- 
gen line with constant width (since we assume local depo- 
sition of positron energy). This systematic error causes the 
discrepancy between the evolution of line width with epoch 
observed and our constant line width in the modelling, and 
it is taken into account in our treatment of temporal line- 
width evolution. Since this temporal evolution is weak, oxy- 
gen and 56 Ni cannot be strongly separated in the observed 
SNe and our modelling approach appears to resemble the 
physical situation quite well (e.g., a very central 56 Ni dis- 
tribution would cause a rapid decrease of line width with 
time). Thus, we are convinced that our description of the 
56 Ni core of the SNe is sufficiently accurate. 

To enable a direct comparison to a parameter which 
can be obtained without any detailed modelling, we calcu- 
lated the half width at half-maximum intensity of the oxygen 
doublet [O i] AA6300, 6364. These two lines are separated 
by ~64 A, which translates to a velocity of ~3000 km s -1 . 
With decreasing density the intensity ratio of the two fines 
will shift from 1:1 to 1:3 ( |Li fc McCray|1992| piugai|1992[ ), 
and in the nebular phase the ratio should lie somewhere be- 
tween 1:2 and 1:3. The error in our characteristic velocity 
caused by the superposition of these two lines will therefore 
be small. 

In Section 2 we already mentioned the advantages of 
both v a and V50. While v a is time consuming to compute, 
W50 is less exact in characterising the velocity of the central 



ejecta. In Figure[3]one can see a rather linear trend between 
both characteristic velocities, as expected. There are, how- 
ever, some outliers. The SNe in the upper-left corner (where 
v a - 1500 km s" 1 >> U50) are SNe 1998bw, 2004gq, and 
20071. Both BL-SNe (SN 1998bw, 20071) have very convex 
line profiles with a broad base and a sharp peak. SN 2004gq, 
the most extreme outlier, also has a convex shape; in addi- 
tion, it shows a kink in the line profile, which explains the 
large difference between v a and vso for this case. The SNe 
in the lower-right region (where v a — 1500 km s _1 << 1150) 
are SNe 1990aa, 1990B, 2005bf, and 2006T. SN 1990aa has 
a prominent "spike" exactly at half height in the blue wing 
of the line profile, artificially increasing U50 without affect- 
ing v a - SN 1990B has a very steep blue wing. SNe 2005bf 
and 2006T both have double-peaked [O 1] profiles with a 
very steeply falling blue wing. The width at half height is 
almost the same as at the base of the line, which explains the 
low v a /vso ratio (see Figure [9] for the different line profiles). 
These examples demonstrate how U50 can give a mislead- 
ing picture of the characteristic core velocity in some cases. 
Hence, while in general ^50 seems to be a good proxy of 
core velocity, it should not be used if the [O 1] profile shows 
broad double peaks, spikes, and kinks, or unusually steep or 
convex wings. 

Both estimates of characteristic velocity are affected by 
the presence of an underlying continuum, which can either 
be the host galaxy or residual continuum emission form the 
SN. It is often difficult to distinguish the continuum from the 
SN spectrum. We tried to overcome this problem setting a 
characteristic minima around the oxygen line to zero flux by 
removing some linear function from the spectrum. However, 
the continuum in most cases is probably not represented by 
a linear function. Thus, some additional flux is almost al- 
ways present in the region of an emission line, causing an 
error in the modelling procedure and in the determination of 
the half height. Consequently, it is not possible to obtain a 
single "best" result for a given spectrum. Depending on the 
quality and the specific shape of a spectrum, there might 
be a variety of possible background subtractions. To cope 
with this problem, we tried to model the extrema of what 
seemed plausible subtractions — of course a rather arbitrary 
approach. We modelled several different SNe in this manner 
to get a quantitative estimate of the typical uncertainty and 
found that an error of ±7% should cover the plausible range 
(e.g., see Figures [l] and [2J where the differences in charac- 
teristic velocity for the models shown are about 5%). 

Another general uncertainty which we cannot quan- 
tify is introduced by possible global asphericities of the SN 
ejecta. In such a case the projected velocities might be con- 
siderably lower than the actual ejecta velocities, so the SN 
kinetic energy may be underestimated. As long as the ejecta 
geometry and inclination are unknown, this problem could 
not be removed by three-dimensional (3D) modelling either. 

For 19 SNe (~35% of our sample), the shell modelling 
approach is not adequate to fit the central parts of the oxy- 
gen doublet. This suggests that at least 35% of the CC- 
SNe of our sample might be aspherical in the very centre 
(for other explanations, see e.g. [Milisavljcv ic et al.||2009"| ). 



Taubenberger et al. ( 2009 1 came to a similar conclusion. Of 



course, asphericities cannot be ruled out for the rest of our 
sample, even if the shell modelling approach was sufficient 
to obtain a good fit to the full line profile. As v a is dom- 
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in this paper). The listing of line widths at different epochs 
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Matheson et al. (20011 seems to confirm that line width 
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Figure 8. The characteristic velocity v a [km s 1 ] against the 
ratio (-E cj ,tot/Af c j,tot) 1/2 [(10 51 ergs/Me) 1 / 2 ], which is a proxy 
for the outer ejecta velocity. There seems to be a weak correlation; 
however, the scatter is large and there is almost no predictive 
power in that relation. 



inated by the outer parts of the line profile, a discrepancy 
between the model and the observation in the central parts 
of the line does not cause large errors in v a . 



5 RESULTS 

Normal SNe of different types seem to have quite similar 
average core velocities v a : 4402 km s _1 (a v = 403 km s _1 ) 
for SNe lib (5 objects), 4844 km s _1 (a v = 935 km s _1 ) 
for SNe lb (13), and 5126 km s" 1 (u„ = 816 km s" 1 ) for 
SNe Ic (27). SNe Ic have an average v a only slightly higher 
(5%) than that of SNe lb and ~15% higher than that of 
SNe lib. SNe BL-Ic have on average a higher v a : 5685 km 
s _1 , and a similar scatter (a v = 824 km s _1 ) to type lb and 
Ic SNe. The sample of SNe BL-Ic comprises 12 SNe — 10 of 
Type Ic, 1 of Type lb, and 1 of Type lib. Our uncertainties 
are rather large, but as long as there is no systematic over- 
or underestimate for one of these groups the ratio of their 
averages should be a reliable quantity. 

In Section 4 we have shown that if there is a physical 
correlation between 56 Ni mass and characteristic core veloc- 
ity it is weak. There also seems to be no correlation between 
total ejecta mass or kinetic energy and core velocity. There 
is only a weak correlation between the ratio of total kinetic 
energy to total ejecta mass and core velocity. On average, 
SNe with higher outer ejecta velocities have higher core ve- 
locities, but this trend is weak and shows a large scatter (see 
Fig ure |8]) . 



Matheson et al. (20011 found that SNe Ic have signifi- 



cantly higher kinetic energy to mass ratios than SNe lb, as 
measured from the line width at half maximum. First we 
compare the SNe contained in Matheson et al. ( 2001[) and 
our sample. We find that the velocity estimates of |Mathe-| 
|son et ah] ( |2001[ > and v^o in this paper agree rather well (to 
within ~10%) for most SNe, consistent with our estimate of 
the general uncertainty. The difference is probably caused 
by the different sizes of the two samples (2 SNe lb and 12 



SNe Ic in Matheson et al. (2001 1; 13 SNe lb and 27 SNe Ic 



is rather constant after the nebular phase has been reached 
(as discussed above). 

Taubenberger et al.| ( p009] | give the FWHM for 39 CC- 
SNe, obtained by fitting the [O i] AA6300, 6364 lines with 
Gaussians. The absolute values and temporal evolution of 
these velocities are consistent with our V50 estimates: almost 
all velocities obtained by Gaussian fitting agree with V50 to 
10% or better, which is within the estimated errors for v$o- 
Velocities obtained with Gaussian fitting are often slightly 
lower than vso, probably because of the contribution of [O 1] 
A6364 to H50. The average velocity for SNe BL-Ic (6 objects) 
from Gaussian fitting is ~ 3670 ± 860 km s" 1 , while the 
average U50 is ~ 4010±600 km s" 1 for these six objects. Since 
we expect U50 to be slightly higher than velocities obtained 
from Gaussian fitting, these estimates are consistent. 

The velocities of SNe BL-Ic are on average only 10% 
higher than those of regular SNe Ic. Interestingly, four BL- 
SNe have core velocities typical of regular SNe. For two of 
them (SNe 1997ef, 2003dh) our spectra are extremely noisy, 
so the errors could in principle be larger than the estimated 
14% (although we do not expect that), which would allow 
higher velocities. On the other hand, for the two other SNe 
BL-Ic with low core velocities (SNe 1997dq, 2003bg) the 
spectra are rather good, and high core velocities can clearly 
be ruled out. Since SNe 1997ef and 1997dq are very similar 
( Mazzali et al.|[2004[ ) in light curve and spectral behaviour, 
it seems likely that our estimate for SN 1997ef (which agrees 
very well with that for SN 1997dq) is also correct. SN 2003dh 
has a very small v a . 

For the GRB-SN 2003dh we might observe a projec- 
tion effect. Given that a GRB was detected, we probably 
observed the SN close to the poles. If the central region 
contained an oxygen-rich disc expanding preferentially near 
the equatorial plane, the projected velocity of the [O 1] line 
would be much smaller for a pole-on view. For example, if 
the disc had an opening angle of 45° , the projected velocity 
would be ~ 70% of the radial ejecta velocity, which would 
mean that SN 2003dh would have been observed to have a 
"regular" (relatively high) BL-SN core velocity had it been 
observed close to the equator. It is not clear whether the 
other three BL-SNe with low v a can be explained in a sim- 
ilar way. One would expect roughly one third of all BL-SNe 
to be observed close to the poles and two thirds close to 
the equator, which is roughly consistent with 4 slow, 7 in- 
termediate, and 1 fast BL-SN (Figure 9 shows that there is 
a rather steep transition from slow to fast between 40° and 
80°). 

As the asphericity can be much weaker in the outer lay- 
ers of the SN than in the centre, this geometrical interpreta- 
tion is not in conflict with the early-time classification. For 
SN 1998bw, |Tanaka et al] ( |2007[ ) have shown that the clas- 
sification as BL-SN (from the early spectra) is not affected 
much by possible ejecta asphericities. 

To test whether this scenario is plausible, we computed 
3D nebular spectra of equatorial oxygen discs with different 
opening angles for several observer inclinations. We modified 
the shell nebular code to operate in 3D. For our simulations 
we use a spherical grid with ~ 6000 cells. We simulate discs 
of 56 Ni and oxygen reaching out to 10,000 km s _1 with open- 
ing angles of 20, 40, and 60 degrees. Varying the observer 
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Figure 9. Left: SN 1990B, whose steep blue wing causes a small 
ratio of v a and vgo. Middle: SN 1983N, which shows an average 
ratio of v a and v^o- Right: SN 2004gq, whose convex shape causes 
a large ratio of v a to 1)50, which is additionally increased by the 
small notch in the blue wing. 





1 










9 







8 


x 





7 



E 





6 


> 





5 







4 


> 











3 







2 







1 











10 20 30 40 50 60 70 80 90 
Observer Inclination 

Figure 10. Simulated evolution of normalised v$o with observer 
inclination (in °, 0° polar, 90° equatorial view) for a 56 Ni/0 
disc with opening angles of 20° (blue, lower line), 40° (green, 
central line), and 60° (red, upper line). Opening angles of only 
20° would cause very low projected velocities, which seem not 
to be observed. Opening angles between 40° and 60° seem to be 
able to explain the observed range of projected velocities, and 
are consistent with the distribution of high and low characteristic 
BL-SNe (see Figure 3). 



angle in 10° steps from the pole to the equator, we calculate 
«50 for all three disc opening angles. As can be seen in Figure 
|10[ oxygen discs with opening angels between 40° and 60° 
would be able to explain the observed difference of BL-SNe 
core velocities and would also roughly agree with the ratio 
of low and high v a BL-SNe cores. Discs with larger opening 
angles would result in smaller variations of the characteristic 
velocity, while very thin discs would cause larger differences. 
Of course, there might be some different explanations (e.g., 
different density profiles of the SN progenitors might cause 
different ratios of outer and inner velocities), but the geo- 
metric explanation is the most straightforward, especially as 
asphericities are expected for a substantial fraction of CC- 
SNe (e.g., |Maeda et al.|2008 |. 

It is important to note that GRB-SNe 1998bw and 
2006aj both have rather high v a . In the scenario described 
above, we argued that GRB-SNe might show low inner ejecta 
velocities, as we would observe them close from their poles. 
For GRB/XRF 060218 there is ongoing discussing whether 
it is a low-energy GRB or an XRF. In the latter case the ob- 
served X-ray emission would allow no conclusions about the 
observer's inclination. Alternatively, and this refers to SN 
1998bw as well, SN-GRBs might have larger opening angles 
than estimated for high-redshift GRBs (the determination 
of GRB opening angles is highly uncertain anyway, since it is 
not clear whether the only known method to do so, by mea- 
suring jet breaks, is reliable). Therefore, these two objects 
might challenge the purely geometrical interpretation pre- 
sented here. Studying the inner ejecta of future GRB-SNe 
might shed some light on this issue. 

Finally, we mention another interesting object in our 
sample, SN 2007bi. Based on our estimate of the 56 Ni mass 
we could speculate that it is a pair-instability SN. The neb- 
ular spectrum of this SN is very different from that of the 
other CC-SNe of our sample. The low core velocity is consis- 
tent with theoretical predictions that PI SNe produce mas- 
sive ejecta at moderate velocities ( Scannapieco et al.|2005 |. 



6 SUMMARY AND CONCLUSIONS 

We estimated 56 Ni masses for 29 SNe for which 56 Ni masses 
were not previously known. Their average agrees with the 
average of 18 CC-SN 56 Ni masses estimated before rather 
well, however there might be a small systematic over- 
estimate of the 56 Ni mass in this work. Individual estimates 
might be quite erroneous given the poor quality of the data 
set. We then measured characteristic velocities for 56 CC-SN 
cores, which range from 3000 km s~ x to 7000 km s -1 (w Q ). 
Several BL-SNe with high-velocity outer ejecta have high- 
velocity cores as well, but some BL-SNe do not. We have 
shown that this might be due to ejecta asphericities, which 
are expected theoretically and might have been detected by 
different methods before. We found that the average core ve- 
locities of SNe lib are slightly lower than core velocities of 
SNe lb and SNe Ic. SNe lb and Ic have very similar average 
core velocities. SNe lib show much less variance (a v « 400 
km s _1 ) of their core velocities than SNe lb and Ic (a v « 850 
km s _1 ). 

There seems to be no strong dependence of core veloc- 
ity on 56 Ni mass. We also checked for correlations between 
total ejecta mass or kinetic energy and the core velocity and 
found none. There is only a weak correlation between the 
ratio of total kinetic energy to total mass and the core ve- 
locity (BL-SNe have the highest core velocities on average, 
but the scatter is so large that there is almost no predictive 
power). Therefore, although the total mass of a SN might be 
estimated (if the spectrum is flux calibrated), it is not pos- 
sible to estimate the SN total kinetic energy from a nebular 
spectrum with good accuracy, since the core velocity seems 
to correlate with the outer ejecta velocities only weakly. We 
will further study the relation between outer and inner ejecta 
velocities in future work. 

The uncertainties in our estimates are rather large. 
They are caused by the general properties of our method 
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(background subtraction, reddening, time evolution of line 
width, uncertainty on epoch, 56 Ni distribution, and 56 Ni 
mass) and are difficult to improve upon when only limited 
high-quality data are available that cover different CC-SNe 
at several different epochs. Therefore more and better data 
is needed. 
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APPENDIX A 

In CC-SNe nebular spectra, the [O I] AA6300, 6364 doublet 
often has a double or even triple peaked emission profile. Re- 



cently, Milisavljevic et al. ( 2009 1 studied these peaks for 20 
CC-SNe and found that they are often separated by ~64 A. 
Since this coincides with the doublet separation, |Milisavlje-| 
vie et al. ( 2009 1 concluded that the multi-peaked line pro- 



files might be caused by some absorption processes, rather 
than by ejecta geometry. In our sample there are 24 SNe 
with clear multi-peaked [O I] AA 6300, 6364 profiles, and we 
measured their peak separation A (see Table 4). 

For SNe 2004ao, 2005kl, 2006T, and 2008ax our results 



agree with those of Milisavljevic et al. ( 2009 1. For SN 2003jd, 



Milisavljevic et al. (2009) did not measure the separation 
of the two largest peaks, while we do, which explains the 
different value given in this paper. 

Although about 50% of the SNe might be consistent 
with a separation of 64 A(within the uncertainties), the 
other 50% are not. For SNe 1985F and 2002ap it seems 
likely that the second (very weak) peak is caused by the 
[O i] A6364 line, since the first peak is roughly at 6300 A 
while the second is at ~6360 A. For SNe 1990U, 1991aj, 
2000ew, 2003jd, 2004ao, 2004gt, 2006T, 20071, 2007bi, and 
2008ax, the two peaks are of similar strength and are cen- 
tered around 6300 A. The separation observed in SNe 1990U, 
1991aj, 2003jd, 20071, and 2007bi is clearly not consistent 
with 64 A. For the remaining 12 SNe the situation seems 



even less clear and we refer to Milisavljevic et al. ( 2009 1 for 
some possible explanations. 

Interpreting the separation of these peaks as a geomet- 
rical effect, one obtains typical velocities between 2000 km 
s _1 and 4000 km s" 1 (1000 km s _1 to 2000 km s _1 when 
considering half width) which is of the same order as the 
characteristic core velocities measured in this paper. There- 
fore the observed clustering around ~ 3000 km s _1 (close 
to the 64 A doublet line separation) might be a coincidence 
caused by the typical velocity of the SNe cores. 

In conclusion, one can say that ejecta geometry remains 
an interesting explanation for split-top line profiles. 



SN 


A [A] (this work) 


1985F 


58 ± 4 


1990B 


65 ± 10 / 65 ± 10 


1990U 


48 ± 10 


1990aa 


45 ± 10 / 61 ± 6 


1990aj 


49 ± 6 


1996N 


65 ± 10 


1996aq 


63 ± 6 


2000ew 


59 ± 6 


2002ap 


58 ± 4 


2003bg 


38 ± 8 


2003jd 


100± 20 


2004ao 


64 ± 10 


2004dk 


72 ± 10 


2004gt 


64 ± 8 


2005N 


50 ± 10 


2005bf 


50 ± 20 


2005kl 


65 ± 10 


20061d 


42 ± 6 / 61 ± 8 


2006T 


70 ± 20 


2007C 


52 ± 10 / 50 ± 10 


20071 


45 ± 10 


2007bi 


47 ± 6 


2008D 


45 ± 10 


2008ax 


61 ± 4 



A [A] Milisavljevic et al. ( 2009 i 



64 ± 5 

65 ± 3 



65 ± 4 

63 ± 3 

64 ± 1 



Table 4. Peak separation measured in this work and by |Mil-| 
lsavljcvic et al.| | [2u09l l. Two values are given if three peaks are 
observed; the separation is then measured from one peak to the 
next. For all SNe but SN 2003jd the estimates agree rather well. 
For SN 2003jd the difference is caused by the fact that |Milisavlje-| 
|vic et ah] ( |2009| ) did not measure the separation of the two largest 
peaks. 
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